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We present a unified approach to the Thermodynamic Bethe Ansatz (TBA) 
for magnetic chains and field theories that includes the finite size (and zero tem- 
perature) calculations for lattice BA models. In all cases, the free energy follows 
by quadratures from the solution of a single non-linear integral equation (NLIE). 
[A system of NLIE appears for nested BA]. 
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We derive the NLIE for: a) the six- vertex model with twisted boundary condi- 
tions; b) the XXZ chain in an external magnetic field h z and c) the sine-Gordon- 
massive Thirring model (sG-mT) in a periodic box of size (3 = 1/T using the 
light-cone approach. This NLIE is solved by iteration in one regime (high T in 
the XXZ chain and low T in the sG-mT model). In the opposite (conformal) 
regime, the leading behaviors are obtained in closed form. Higher corrections can 
be derived from the Riemann-Hilbert form of the NLIE that we present. 



1. Introduction 

The computation of thermodynamical functions for integrable models started 
with the seminal works of C.N. Yang and C.P. Yang [1], of M. Takahashi [2] and 
M. Gaudin [3]. In Refs. [2,3] the free energy of the Heisenberg spin chain is written 
in terms of the solution of an infinite set of coupled nonlinear integral equations, 
derived on the basis of the so-called "string hypothesis" . 

These equation are commonly known as Thermodynamic Bethe Ansatz (TBA) 
equations. Nowadays, the TBA is mostly used as a nonperturbative tool to inves- 
tigate Perturbed Conformal Field Theories and has been considerably generalized 
(see e.g. [4,5]). In all cases some some kind of string hypothesis is used. 

In ref.[6], we proposed a simpler way to solve the thermodynamics of the XXZ 
chain and of the sine-Gordon field theory by means of a single, rigorously derived, 
nonlinear integral equation. This method can be applied to a wide class of models 
solvable by Bethe Ansatz, and can provide also the excited-states scaling functions 
of the sG field theory [7]. Similar methods were used to also determine some 
correlation lengths at finite temperature in the XXZ chain [8]. 

The purpose of the present paper is threefold. First we generalize the results 
of ref.[6] in the presence of external magnetic field for the XXZ Heisenberg model. 



Second, we show that this approach provides an unified method to compute 
both at finite temperature and infinite size and at zero temperature and finite 
size. Third, we recast our NLIE in two alternative forms: a) as an infinite set of 
algebraic equations of BA type and b) as a Riemann-Hilbert problem. These forms 
are appropriate to compute the asymptotic behaviors around the conformal regime 
(low T for the magnetic chains and high T for field theory). In the recursive regime 
(the opposite to the conformal one), our NLIE can be easily solved by iteration. 

Continuum quantum field theories at finite temperature T are obtained from 
periodic boundary conditions (PBC) in the imaginary time with period (3 = T _1 . 
By Wick rotation, this leads to a periodic spatial box of finite size f3. What we 
achieve here is the analogous to a Wick rotation on the exact Bethe Ansatz (BA) 
solution on the lattice. 

Our starting point is the BA diagonalization of a row-to-row inhomogeneous 
and twisted transfer matrix t(X, #i,#2> • • • ,9n,uj) where there is an inhomogene- 
ity 9 n attached to each site and a twist u in the boundary conditions (ui = 
corresponds to PBC). 

Setting the inhomogeneities 9 n = yields the homogeneous six-vertex transfer 
matrix. 

Choosing the inhomogeneities alternating 9 n = (—l) n ~ 1 9 yields a diagonal-to- 
diagonal transfer matrix Tjy(9) = t2iv(6 l , 9, 0) [20]. For 9 — ► 0, T^{9) gives at order 
9 the XXZ hamiltonian. 



Tn{9) 6 =° 1 - .ry—Hxxz + 0(9 2 ) (1.1) 
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and hence, 
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This relation permits to write the XXZ partition function 

Z = Ti f e -P Hxxz \ 

in terms of the diagonal-to-diagonal transfer matrix Tl(8). Moreover, using the 
crossing symmetry of the statistical weights one can prove that Z in the infinite 
volume limit is given just by the ground state eigenvalue of Tjy(9) in the iV — ► oo 
limit. This alternative way is much simpler than the usual TBA approach, where 
one has to sum over an infinite number of states with an arbitrary number of holes 
and an arbitrary number of strings of any length. What we do here is analog to 
duality in S-matrix theory, where the sum over an infinite number of resonances 
in the s-channel can be replaced by the exchange of one reggeon in the t-channel. 

Setting 9 = —iQ with real, we get the lattice field theory evolution operators 

in the light-cone approach [13,20] 

U R (Q) =a(-2i&)- N t 2N (-i&,-i&,0) , U L (Q) = a(2iQ) N t 2N (+iQ, -i®, 0) _1 

where a(±2iO) is a unitarization factor. That is, for imaginary 9, we interpret 
the lattice as a discretized Minkowski spacetime [15] and hence the operators 
Ur and Ul describe space and time translations in the light-cone directions of a 
discretized Minkowski spacetime. If we call S the lattice spacing, 

U = e~ iSH , U R XJ L X = e iSP (1.2) 

define the lattice hamiltonian (H) and momentum (P) operators. 

Moreover, the operator i2iv(A, — iQ, 0) is the generating functional of local 
conserved charges (including H and P) for |ImA| < 7/2 [21]. In the region 
— (it — 7/2) < ImA < —7/2, i2iv(A, — iQ, 0) generates non-local conserved charges 
[21]. 



We diagonalize the transfer matrix t(X, 9\, 02, ■ ■ ■ , @N,w) by algebraic Bethe 
Ansatz (BA). The BA the expresses eigenvectors and eigenvalues in terms of the 
roots of M coupled algebraic equations (1 < M < N), the so called Bethe Ansatz 
equations (BAE). 

Using contour integral techniques, we rewrite these BAE as a non-linear inte- 
gral equation. This equation may take different forms. A compact form is 

+00 
Z{\i) = y?(/i) + 2 Im / d fi' G(pt - »' - irj) log [l + e iZ ^' +irl) ] (1.3) 

—00 

where the kernel G(X) is explicitly given by eq.(2.24) and 
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(1.4) 

Z([j,) is the unknown in eq.(1.3). The real parameter 77 can be freely chosen in the 
interval < rj < 7/2 . 

We have unified three different problems (finite size vertex model with TBC, 
the XXZ thermodynamics in a magnetic field and the sine-Gordon field theory in 
a finite volume) into the single NLIE (1.3)-(1.4). 

Once we have the solution Z(\) of eq.(1.3), we can compute the transfer matrix 
eigenvalues by quadratures. We consider three relevant applications 

a) The finite size (homogeneous) six-vertex model with twisted boundary con- 
ditions (TBC). 

b) The (infinite) XXZ chain at temperature T = /5 _1 in an external magnetic 
field h. 

c) The massive Thirring-sine Gordon model (mT-sG) on a periodic box of size 

(3. 



In the first case, we find for the six vertex partition function on a double 
periodic square lattice (T x 2N ) 



lim T logZ, 
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-2Nfo(\,u) + L 2N (\,uj) 



1.5) 



For large N, we obtain that Ln follows the conformal behaviour for TBC with 
unit central charge [see sec. 7.1]. Namely, 



T N^oo 7T ( T\\ 

L2N = -^t-j tan — 
bi\ \ 7 
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In the case b) , we relate the XXZ free energy at temperature T and on an external 
magnetic field 2h in the ^-direction with the maximal eigenvalue of the alternating 
transfer matrix £2jv(±#, &, —i/3h) with TBC. Notice that the twist, —iflh, is purely 
imaginary for real magnetic fields h. This new result allows us to express the XXZ 
free energy f(/3, h) in terms of the solution Z(X) of our non-linear integral equation 
(NLIE) (1.3). 



f(/3,h) = /i + 2Jcos7 + 



sm7 



d\ 



log 



1 + e ~ iZ ^ 



2iri sinh A sinh(A — -ry) 



1.6) 



Finally, we use the light-cone approach [13, 20, 15, 21] to deal with c). We write 
the energy in the mT-sG model on a ring of length (3 in terms of the alternating 
six vertex transfer matrix [£(=pO, —iQ, 0)] eigenvalues. The basic relation is the 
the formula for the unit evolution operator for lattice field models, 
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where 5 is the lattice spacing, plays the role of UV cutoff and a(±2«0) are 
unitarizing c- number factors. In the continuum limit, S — > 0, N — > 00 with 



(3 = NS fixed. At the same time tends to infinite as 
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We find for the scaling function 



E{{3) = Jim [E N - E c 
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where Z(X) satisfies [cfr. eq.(1.3)] 
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Z(X) = m/9sinh(7rA/7) + 2 Im / dfiG(X — /j, — irj) log 
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provided we choose < r) < 2min(7, n — 7). 
Shifting A by 17/2, eq.(l.lO) takes the form 



ey = m/3 cosh 9 — Gq* Lf + G\* Lj 



XIV, 



where 

e/(0) = -iZ(j0 + ry/2), e f -=ej 

L f = log(l + e e , L 7 = log(l + e £ /) = L~ f (1.12) 

G o (0) = £G( J0) , Gi(0) = G o (0 + in - zO) 

The kernels Gq(6) and Gi(#) are just the 6? derivative of the fermion-fermion and 
fermion-antifermion physical phase-shifts. This transparent physical interpretation 
of eq. (1.11) suggests the form that it should take in other field theory models. 

Now we face the problem of solving NLIE like (1.3) or (1.10) in order to get 
physical results. It must be noticed that we find just one NLIE for finite temper- 
ature and finite size situations. This situation contrasts with the standard TBA 



which for nonrational values of 7/7? involves an infinite number of nonlinear inte- 
gral equations for the infinitely many different types of "magnons" describing the 
energy degeneracies of the physical fermions and antifermions [23] . Therefore, our 
equation effectively provides a resummation of the magnon degrees of 
freedom. 

There are essentially two regimes in our NLIE. We call them 'recursive' and 
'conformal'. In the 'recursive' regime the NLIE can be simply solved by iterating 
the inhomogeneity. This happens for small j3 in the XXZ chain and for large j3 in 
field theory. 

The 'conformal' regime is the opposite to the 'recursive' regime. In it, the 
calculations are more subtle. The dominant behaviour (large /3 for the XXZ chain 
and small /3 in field theory) , follows in closed form with the help of the Lemma in 
sec. 7. This lemma avoids the detailed resolution of the NLIE for the dominant 
contributions. Higher corrections follow using the Riemann-Hilbert form of our 
NLIE (sec. 7.4). We find that the energy in the sG-mT model has a series of 
terms with (m/3) -1 times integer powers of (m/?) 47 '^ plus the central charge term 
proportional to (ra/3) -1 , in agreement with the predictions of Perturbed Conformal 
Field Theory. 

Since the inhomogeneity in the NLIE for field theory is m/3 cosh 8, it is useful 
to define for small m(3 kink solutions, e&(0) of the NLIE with inhomogeneity e 
[see eq.(7.19)]. We obtain in this way for the energy E(f3) at high temperatures 



E P =° ~ - ^ cot [7r 2 /2 7 ] - £ J dd coshtf [£(9) + £(9)} (1.13) 



The first term is the universal conformal Casimir energy, from which we read the 
correct central charge c = 1 of the mT-sG model. The second term, linear in /3, 
exactly coincides with minus the scaling bulk free energy [see eq.(5.7)]. The last 
integral in eq.(1.13) represents the resummation of the perturbed conformal field 



theory around the massless Thirring model. The fact that it must contain non- 
integral powers of r, causing non-analyticity at r = 0, can be established directly 
from the original equation (7.18), since this cannot be expanded in a Taylor series 
of r. 

In the 'recursive' regime, we derive in sec. 6.1, the high temperature expansion 
for the XXZ chain in the external field h. One only needs the Cauchy theorem 

to evaluate the high T expansion coefficients. The first five terms are explicitly 
computed in sec. 6.1: 

/(/?) = - p" 1 In 2 + J cos 7 - (3 [J 2 (l + ±cos 2 7) + \h 2 ] + 1 cos 7 ( J 3 + h 2 ) 



- (3 S [(\ + icos 2 7 - ^cos 4 7 )^ 4 - J 2 h 
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which indeed agrees with the high T expansion. 

In sec. 6.2 we obtain the ground state energy for the mT-sG model for small 
T. It is non-analytic at T = 0. We find up to contributions 0(e~ 3m ^) 
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(1.14) 

where K\ ix) stands for a modified Bessell function and the constant K^ is defined 
by the integral (6.17). The first term in eq.(1.14) describes free massive bosons, 
the subsequent terms take into account the interactions. 

In sees 4 and 5, we show that our NLIE (1.3) and (1.10) are equivalent to 
an infinite set of algebraic equations of BA type. These equations differ from the 
usual BA equations [14, 15] in two essential aspects: i) they are infinite in number, 
ii) the roots are discretely spaced although they describe models with an infinite 
number of sites. The root spacing tends to zero only in the conformal limit. 

We would like to remark that, contrary to the traditional thermodynamic Bethe 
Ansatz [2,3], our approach does not rely on the string hypothesis on the structure 
of the solutions of the BA equations characteristic of the Heisenberg chain. This 



makes our approach definitely simpler. Most notably, the whole construction of 
the thermodynamics no longer depends on whether 7/7? is a rational or not, unlike 
Takahashi approach [2]. This applys equally well to the problem of the ground 
state scaling function of the sine-Gordon field theory, since the standard TBA 
approach requires the string hypothesis. 

Generalizations of the present TBA approach to higher spin chains as well as to 
magnetic chains and quantum field theories associated to Lie algebras other than 
A\ (nested BA solutions) are relatively straightforward provided crossing symmetry 
holds. 

Our NLIE has been recently used with success in ref.[24]. 

2. Nonlinear Integral Equation 

Let us consider a horizontal line formed by N sites, and associate to each site 
n — 1, 2, . . . , N a triplet of (vertical) Pauli matrices (a^, an, cr%), with standard 
commutation rules 



(T (7 



£<> Vm,n t <J n 



As customary, we also introduce the auxiliary, or horizontal, two-dimensional 
space, labeled by the index n = 0, on which the matrices (<7q, <Tq, Oq) act. For 
any complex value A of spectral parameter and an arbitrary set of inhomogeneities 
0i,$2, ■ ■ ■ ,&Nj the local vertices L n are then written as 

L n = R 0n (\ + O n )Pon (2.1) 

in terms of the spin 1/2 R— matrices of the six vertex model, 

«) — + *^r« + \ KX + °l<) ( 2 - 2 ) 

and of the exchange operators, 

Pkn = o (1 + ^fe ' °n) 



10 



The dependence on 9 of R^ n (9) is contained in the the Boltzmann weights a, b, c, 
for which we take the standard trigonometric parametrization 

a = o (0) = sin(7 - 9) , b = b(9) = sin9 , c = shi7 (2.3) 

where 7 is the anisotropy parameter. Throughout this paper, unless explicitly 
stated otherwise, we shall confine our attention to 7 real and in the interval (0, 7r/2) 
(by periodicity, the complete relevant interval is < 7 < n). 

The inhomogeneous and twisted six-vertex transfer matrix associated to this 
N-sites line reads 

t(X, 0i, 9 2 , ... , 9 N , 00) = e iua ° tr \L X L 2 ...L N ] (2.4) 

The twist angle u defines the relation between spin operators after a translation 
by N sites: 

<+N ± w'n+N = e±lUJ K ± i° V n \ , <+N = < ( 2 -5) 

Notice also that Rkn is t/(l)— invariant, 

e W°i+0 Rkn ( d } = Rkn ( d } e iH* z k +<) (2.6) 



so that the transfer matrix commutes with the z— projection of the total spin 

1 

2 



cz 1 v^A^ z 

D — 2 Z^l a n- 



The diagonalization of t(X, $1, 9 2 , . . . , 0jv, a;) is achieved via Quantum Inverse 
Scattering Method (also called Algebraic Bethe Ansatz) (ref. [14,15]): the eigen- 
states with S z = N/2 - M, M = 0, 1, 2, ... , [N/2], are labelled by the M distinct 
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roots {Ai, A2, • • • , Am} of the Bethe Ansatz Equations (BAE) 



n sinh(A i + i6 n + ij/2) _ ^_ 2iu yr sinh(A j - A r + ij) 
wah(Xj + i$ n -i'y/2) ' AA sinh(Aj - A r - 17) l ' j 



while the corresponding eigenvalues r(A, 61,62, ... , 9n, u) read 



J. J. J. J. cir 



n—1 r 

N M 



^f- sinh(i7/2 + Aj — iA) 



.'2.8) 



+ «-~n^« n rS^;g 



n=l 



We shall primarily be interested in one particular solution with S* z = of the 
BAE, the so-called antiferromagnetic ground state (a.g.s.) solution, which will be 
characterized shortly. Since the request S z = forces N to be even, we replace N 
by 2 A" throughout. Moreover, for our applications, we shall need to consider only 
a special type of inhomogeneity structure, the alternating one 



6 n = {-l) n - l 9/2 (2.9) 

where 6 will either be real, with < 6 < 7/2 , or purely imaginary. Correspond- 
ingly, the object of our interest is the (twisted) alternating transfer matrix 



t 2N (\, 6, u) = t(X, 0i = 6/2, 6 2 = -6/2, ...,6 2N = -6/2, uo) (2.10) 

As is well known, this object is a diagonal to diagonal transfer matrix with 
N links in the horizontal direction [20]. Suppose now that N distinct roots 
{Ai, A2, • • • , Ajv} are given; with them we construct the so-called counting function 

12 



[15] 
Z N {\) = N [0(A + z0/2, 7 /2) + 0(A - i$/2, 7 /2)] - ^ 0(A - A fe , 7) - 2cu (2.1i; 



TV 



fc=l 

where 

0(A,x)^zlog Sin ^ X + ^ , (0(0,*) = 0) (2.12) 

sinn(« — A) 

has the cut structure chosen so that it is analytic in the strip |ImA| < x. The a.g.s. 
solution is now specified by the property 

Z N (Xj) = (-N + 2j - 1)tt , j = l,2,...,N (2.13) 

and enjoys the symmetry (compare with the BAE, eq.(2.7)) 



Xj(lo) = Xj(uj) = -X N - j+1 (-uj) (2.14) 

The existence of this solution is most simply established by numerical calculations, 
for N up to the thousands. The uniqueness will become apparent in the sequel. 
For uj real, all Xj are also real and the a.g.s. can be characterized as the unique 
solution of the BAE with N real roots. If Im uj > (< 0) then Im Xj > (< 0) and 
the a.g.s. goes continuously, by construction, into the real solution as Imw ^ 0. 
Notice that the symmetry (2.14) implies for the counting function itself 



Zn(X,uj) = Zn(X,uj) = —Zn(—X,—uj) (2-15) 

which will play an important role in the applications. Notice also that all these 
last statements really hold true for 7 in the extended interval (0, tt). 

To gain some understanding on the distribution of roots in this a.g.s. solution 
one can try to approximate it in the limit of N large and fixed (i. e. N— independent) 
9 and uj. The idea is the traditional one that the BA root distribution becomes 

13 



the continuous density c c (A) (normalized to 1) in the limit iV — ► oo. Then the 
sum over the roots in the definition of Zjv(A), eq.(2.11), can be approximated by 
N times the integral over the density cr c (A). At the same time the property (2.13) 
implies that 

Jim ^jUiv(A) = 2na c (\) (2.16) 

so that one finds the following linear integral equation for cr c (A) 

+oo 

27ra c (A) = / (A + ^/2,7/2) + / (A-^/2, 7 /2)- f dpi <f>' '(A - ^7) a c (pt) (2.17) 



where <f>'(\) = d<p(\)/d\. Fourier transforming leads to the explicit solution 



cr c (A) = z' c (\) , z c (X) = -tan : 

7T 



sinh(7rA/7) 



(2.18) 



_cos(tt^/[27])_ 

yielding to leading order Zjy(X) ~ 27rNz c (X). Notice that no dependence on the 
twist parameter u is left over in this continuum approximation. One can attempt to 
recover u> by adding to 27rA^ c (A) the constant function /o(A) = — 2u, as suggested 
by the definition itself of Zjy(A), eq.(2.11). However, this continuum approximation 
need not control the 0(1) terms. In fact, to solve (2.17) one has to invert the 
convolution (1 + K) * a c , where K is defined as 



+00 



(A'*/)(A)= / ^(A-x, 7 )/(x) 



(2.19) 



This implies that the correct constant to add is (1 + K) * /o, that is 



2-kNz c (\) - -^— ~ 



7T — 7 

Therefore in the continuum approximation the roots are given by 



A- ~ 

A 3 — 



7 



sinh 



71 



TT0 TX /n . CO , 

cos — cot —77(27 — 1 H 

7 2A^ V J 7T-7' 



(2.20) 



(2.21) 



Although imprecise, this approximation is useful, for numerical tests show that 
the actual values of Aj are very close for arbitrary values of 9 and u. Thus we see 
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that the roots lay on a smooth curve entirely contained in some horizontal strip 
D of the the upper (lower) complex plane for Imcj > (Imw < 0). In turn, T> is 
included inside the strip ImA < 7/2. Moreover, from a host of numerical checks 
as well as from arguments based on the unicity of the a.g.s. configuration, we 
expect to find no other point, in the strip |ImA| < 7/2, where e N takes the value 
(— 1) +1 , besides the roots Xj. We shall now look for the modifications to the 
linear eq.(2.17), which are necessary to go beyond the continuum approximation. 

For sake of simplicity, we shall consider the case N even (after all N — *■ 00 in 
physical applications) and ImcJ > (recall the property (2.15)). By construction, 
the function exp[iZj^(X)} is analytic in the strip (—7 + 6)/2 < Im A < 7/2 and the 
sum in the last term of eq.(2.11) can be written as a contour integral 



TV 



$>(A-A J?7 ) 
3=1 



^0(A-/i,7);flog[l + e^ 
2m a/j, 



'2.22) 



where the closed contour T encircles counterclockwise all the roots Xj (recall that, 
by construction, exp[iZ^(Xj)] = —1). We can take T to run straight from left to 
right at Imz — — rj-, < ??_ < (7 — Re#)/2 and come back straight at Imz — 77+, 
max(ImAj) < i]+ < 7/2 (of course the exact values of rj- and 77+ do not matter 
because of the analyticity and the behaviour at real infinity of e JV ). 



Inserting eq. (2.22) ineq.(2.11), using 1+e N 
by parts yields 



JZi, 



-iZ h 



and integrating 



+00 



Z N (\) + 



dx 

2^ 



<f>'(\ - x, i)Z N {x) = N [0(A + 16, 7/2) + 0(A - iB, 7/2)] - 2u 



-1 I 7T 4>{^ ~ x ~ iV+, 7) log 

2n 



1 + e 



iZ N (x+irj + ) 



+00 

/dx 
— 4>'(X - x + ??7_,7)lo^ 
2it 



I _|_ e -iZ N (x-ir)-) 



(2.23) 
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Notice that the function exp[— iZjq{X)\ is analytic in the strip (7 — Re#)/2 > 
ImA > —7/2, so that ??_ may be chosen to vary up to 7/2. We next introduce the 
convolution operation 



+00 
1 f dk sinh(7r/2 - 7) A; ;/iA 



G = K*(1 + Ky i , G(X)= / t^— T7 Vt-tt; — , , , ^ e ( 2 - 24 ) 

f 47rsinh(vr-7)A;/2 cosh(7/c/2) v ' 



so that eq.(2.23) can be rewritten 



+00 

1 _|_ e iZ N (x+iri + ) 



Z N {\) = 2nNz c (X) — i dx G(X -x- ir} + ) lo 

ix -7 J 

—00 

+CX) 

-iZn(x— if]-) 



+ ij dxG (X- X + i V .)lo, 



l + e 



(2.25) 
This expression is exact and exhibits (although only implicitly) all the corrections 
to the continuum approximation (2.20). If Imw = 0, so that the roots Xj are real, 
the property Zn(\) = Zjv(A) allows to interpret the identity (2.25) as a nonlinear 
integral equation for the counting function itself. In fact, setting r)+ — r\- = 77, 
Zjv(A + irj) = ie(X), G(X + iy) = G y (X) and 

g(X) = -2niNz c (X + irj) + - — . (2.26) 

1 — 7/7T 

we obtain, in compact notation 

e = g - G * log(l + e- e ) + G 2v * log(l + e~ l ) (2.27) 

We recall that 77 can be chosen freely in the interval < rj < 7/2, since the original 
assumption < 7 < n/2 implies that G(X) is analytic in the strip |ImA| < 7. Eq. 
(2.27) may nontheless be extended to the other region n/2 < 7 < n, where G(X) 
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is analytic for |ImA| < n — 7, simply by requiring that < 77 < g-min^, n — 7). 
For 77 — > 7/2, eq.(2.27) agrees with the nonlinear integral equation derived in ref. 
[18] by different methods. 

To convert eq.(2.25) into a constructive integral equation also when Imou 7^ 



we can exploit the symmetry (2.15), namely Zn(\, u>) = Zn(\, uj) and combine the 
two cases Imw > and Imw < 0. In fact, the change uj —^ uj simply amounts to 
change the restrictions on r]±, which now read < 77+ < 7/2 and — max(ImAj) < 
i]- < 7/2. It is again convenient to fix 77+ = r\- — 77, although the lower bound 
for rj actually depends on the roots and introduces an extra nonlinearity in the 
problem. The choice 77 = 7/2, however, (recall that we assumed < 7 < 7r/2) 
would eliminate this subtlety. Then we set 



Zn(\ + irj; Re a; ± ilmw) = ie±(\) 



a 



2ttNz c (X + it]) (Reu ± ilmw) = ig±(X) 

71 — 7 



(2.28) 



and arrive at the system of two coupled nonlinear integral equations 

e± = g± - G * log(l + e- £± ) + G 2v * log(l + e~^) (2.29) 

Let us recall that for lmw = 0we can choose 77 — > ^min(7, tt — 7), while for nonzero 
Imu and 7 < n/2 we can fix rj — ► 7/2. This choices are assumed from now on, 
unless explicitly stated otherwise. 
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3. Finite size corrections 

The nonlinear integral equations (2.29) are exact for any N. In the large N 
limit they allow to find the finite size corrections to the eigenvalue (2.8) of the 
transfer matrix for the a.g.s. configuration. For simplicity we shall consider here 
the homogeneous case (i.e. with 6 = 0) which is relevant in the calculation of 
the partition function Zqv of the standard six-vertex model on a very long torus. 
Indeed on a doubly periodic square lattice with T sites in the vertical direction 
and 2N sites in the horizontal direction Zqv reads 



Z 6V = Tr[t 2N (X,0,u)] T (3.1) 



Actually, this periodic lattice has a vertical seam across which the vertical statis- 
tical variables a = ± acquire the factor e , as evident from eq.(2.5). To simplify 
the discussion, in this section we shall take uo to be real. 

In the cylinder limit T — > oo at fixed N, the largest eigenvalue T max of the 
transfer matrix dominates and we have 



lim T 1 logZ 6V = logT max = -2Nfo(\,u)+L2 N (\,uj) (3.2) 

T— >oo 



where /o is the free energy per site in the infinite N limit, and L 2 n represents 
the finite size corrections. With our original choice of a real anisotropy 7, the six 
vertex model is in the critical regime, so that we expect L 2 n to be governed for 
large N by conformal invariance. 

It is well known that the T max exactly corresponds to the a.g.s. configuration 
analyzed in the previous section, so our problem is now to relate T max to the 
solution of the integral equation (2.27). To this end we first notice that, provided 
< Re A < 7/2, the second term in the expression (2.8) for the a.g.s. eigenvalue is 
exponentially small relative to the first one (for an explicit verification see eq.(6.29) 
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of ref. [21]), and may therefore be dropped. Thus we can write 



log r max = iu + 2N log o(A) + E 



N 



(3.3) 



where 



N 
3=1 



(3.4) 



The sum over the roots can be transformed into integrals by the same procedure 
used in sect. 2 for the counting function, with the result (we must choose here 
r/_ < 7 /2 - Re A) 



+00 +00 

E N = -i /^0 / (/i + ^A, 7 /2)Z iV (A)- fptfiji + iX + iri+n/qioz 



X _|_ e iZ N(v+im 



+00 



+ / 2^^(A* + *A-i77-,7/2)log 



1 + e 



-iZ N (p-ir]-) 



(3.5) 
Eq. (2.25) may now be used to eliminate the term linear in Zjv(A), yielding (we 
choose 77+ = r/_ = i] as usual) 



+00 



En = E c — I dpL < ct c (/jl + iX + irj) lo 



1 + e 



iZ N (fx+irj) 



cr c (/i + A — ir/) log 



1+e 



-iZ N (fi-iri) 



(3.6) 



where _E C is the continuum approximation 



E c = -iN dfi </>(fjL + U, 7/2) o-c(^) 



+00 



:> v ^ dA; sinh[(7r - 7 )fc] sinh(2AA;) < 

fc sinh(7r/c) cosh(A;7) 



(3.7) 
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and (recall eq.(2.18) and that 9 = here) 



aJfi) = -i — = J-[(l - G) * ; ](/i) (3.8) 

vp; 7 cosh(vr ; u/7) 2vr LV ; rjVP ' y ' 

Comparing eqs. (3.2) and (3.6) we obtain 

Zo = -^£ c + loga(A) 

[notice that E c /N is iV-independent] and the finite size corrections expressed in 
terms of the solution of the nonlinear integral equation (2.27) : 

1 + f°, / log[l + e- e M] log[l + ejg] 1 

Lat(A, a;) = 2cu / d/x < — -p- , . . . — r - , , . . > 

7 J I cosh[7r(/i + i\ + w?)/7j cosn|7r(/x + «A — %"i])j^\ I 

— oo ^ ' 

(3.9) 
In section 7.1 we compute the dominant large N behaviour of Ln(\,u) which 
yields the central charge. 

4. Thermodynamics of the XXZ spin chain 
The Hamiltonian of the periodic XXZ spin chain with 2L sites reads 

2L 
HxXZ = ~JJ2 K ff «+1 + a n a l+l - COS 7 {«+! + 1)] (4-1) 

71=1 

and is very simply related to the transfer matrix Tl($) of the symmetric six-vertex 
model on a diagonal lattice with 2L links in the "space" direction. In fact Tl(9), 
with a suitable normalization, can be written 

Tl(0) — C -R2 3-R4 5 •• -R-2L 1-Rl 2-^3 4- • -R2L-1 2L (4-2) 

with the R— matrices Rkn given in eq.(2.2), and for small 9 one finds 

Tl{9) e =° 1 - ttj^—Hxxz + 0{9 2 ) (4.3) 

2 J sin 7 

Notice that the higher order local operators which follow from the 9 expansion 
do not commute among themselves nor with Hxxz, since diagonal-to-diagonal 
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transfer matrices do not commute for different values of 9. This is the reason why 
it is preferable not to call 9 spectral parameter in this context. 

We now observe that eq.(4.3) allows to write an euclidean path-integral relation 
between the classical partition function of the six-vertex model and the quantum 
partition function of the XXZ model [10,11, 12]. Let us consider the XXZ free 
energy per site, /(/3, h), in the presence of an external constant field 2h coupled to 
the conserved z— projection of the total spin S z = \ z~2 n a rv 

f(/3,h) = -\ lim -!-log [lY L-P{Hxxz-2hS z )Y\ ( 44 ) 

p L^oo 2L LI J . 

Setting (3 = pjsin.'y, from eq.(4.3) we then read 

e -/3H xxz = Um Tl (2(3/n) 

Hence we can write 

f(P, h) = ~ lim J- Jim \ogZ LN (2(3/N, h) (4.5) 



where 



Z LN (9,h) = Tr [e^ hSZ [T L {9)\ N ] (4.6) 



is the six-vertex partition function on a square periodic diagonal lattice with 2NL 
sites, with the insertion of an horizontal seam along which each link variable a — ± 
is multiplied by e^ . Notice that any two neighboring elements of this seam can 
be freely moved through a common vertex thanks to the U(l) invariance of the 
R— matrices, eq.(2.6). The two limits in eq.(4.6) cannot be interchanged since the 
degeneracy of T^(0) = 1, that is 2 , is strongly L— dependent. However, the 
numerical value of Zln(9, h) does not change under a rotation by n/2 of the entire 
lattice, nor under the flipping from ± to =p of the value of all link variables on, say, 
the left-oriented diagonals. On the other hand these two operations combined are 
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equivalent to the substitution 9 — ► 7 — 9 in each local i?— matrix i?j & (this is called 
"crossing symmetry", and 7 is then identified as crossing point, see e.g. ref. [9]), so 
that we can write 

Z LN {9, h) = Z NL (l - 9) h = Tr [7M 7 - 9) L h ] (4.7) 

where T^{9)j l is the twisted diagonal-to-diagonal transfer matrix 

TN{Q)h = c *~ e ^-^2 3-^4 5 •• -R2Nl£~ Gl R\ 2-^3 4- • ■ R-2N-12N (4-8) 

Combining eqs. (4.5) and (4.7), we now obtain 

/(/?, fc) = "^ Km -^ Jim log^ L ( 7 - 2^/iV), (4.9) 

It is well known [20] that Tpj(9) is simply related to the row-to-row (untwisted) 
alternating transfer matrix t(X,9,ou = 0) introduced in section 2 

T N {9) = [a(9)a(-9)/c 2 ] N t(9/2, 9/2, 0) t(-9/2, 9/2, 0)" 1 (4.10) 

The presence of the external field can be related to twisted boundary conditions 

[eq.(2.5)] in absence of magnetic field. We consider a purely imaginary twist uo = 
—if3h, and find quite simply 

T N (9) h = [a{9)a{-9)/c 2 ] N t{9/2, 9/2, -i(5h) t{-9/2, 9/2, -i^hy 1 (4.11) 

Therefore the eigenvalues of Tj^{9)y l can be read off the general result (2.8) of the 
Algebraic BA, with the substitutions 

\ = 9/2, 9 n = (-l) n_1 0/2 , uj = -i/3h (4.12) 

Notice that the second term in the r.h.s. of (2.8) then vanishes, so that we have 

, rrr, /m , ,„n u tt sinh(A 7 - + 19/2 + ry/2) sinh(A 7 - %9/2 - i 7 /2) 

eigenvaluejTiv (9) h \ = a(9) N TT — -±-1 4; ^rr — rrr 1 tt, ^rr 

6 i v jnj u 11 sinh ( Ai + ^/ 2 - z 7 /2) sinh(A i - %9/2 + i 7 /2) 

(4.13) 

For < 9 < 7 , the largest eigenvalue A 1 J^ ax (9) of Tj^(9)h is nondegenerate and 

correspond to the solution of the BAE with M = N roots which we have called 
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a.g.s. configuration in sect. 2 . We see therefore that, for h > (h < 0) the roots 
lay in the lower (upper) complex half-plane and are real for h = 0. Moreover, by 
eq.(2.14), Xj(h) = Xj(—h) and Xj = —Xn-j+i- Thanks to the latter symmetry the 
corresponding eigenvalue A^ ax (6)h remains real, while the former will ensure that 
f(/3, h) = f(j3, —h), as required by the spin inversion symmetry of the XXZ chain. 
Finally, the nondegeneracy of hJ^ ax {9)} l implies that the two limits in eq.(4.9) 
commute [11], so that the limit L — *■ oo selects A^ a;r (7 — 2(3 /N)^ and one finds [6] 



f(P> h ) = -™ J im lQ g A N ax (l ~ W/N) h 



2(3 N 



—— lim 

2/3 N^oo 



E N -2Nlog 



sm(2(3/N) 
sin 7 



(4.14) 



where 



jV 



e (A) = i<f>(Xj + i 7 /2, 7 /2 - p/N) - i<j>(\j - z 7 /2, 7 /2 - (3/N) 



(4.15) 



By the same trick used in sec. 2 and 3, we transform the sum over the roots into in- 
tegrals and then use eq.(2.25) to eliminate the term linear in the counting function. 
We obtain in this way 



+ 00 



En = E c — i dX iPn{X + ii]+) 1°£ 



-oo 
+oo 



+ i I dXipN(X — ir)-) lo£ 



1 + e 



1 + e 



iZ N (X+iri^ 



(4.16) 



-iZjv(A— it)-) 



where E c is the continuum approximation 



+oo 



+oo 



E C = N I dXe (X)a c (X) = -2N / — 



dk sinh(7r — 7 )/c sinh( 7 — 2(3 /N) k 



sinli7rA;cosh 7 A; 



(4.17) 
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and 



MX) = -L[(l - G) * e' }(\) 2 sinh ^ A /7) cos(^/ 7 iV) 



2vr LV " 7 cosh(27r A/7) - cos(2tt (3 / 7 N) 

Let us recall that the counting function -Zjv(A) fulfills now eq.(2.25) where, in the 
source term 2nNz c ( y \), we must set 9 = 7/2 — J3/N . We then see that the root 
density of the continuum approximation 

a (X) = 2 cosh ( 7rA /7) sm(ir(3/jN) 

7 cosh(27rA/7) - cos^np/jN) 

is strongly peaked at A = for large N, reflecting the singularity which develops at 
the origin in the N — > 00 limit of the source term in eqs.(2.17) when 9 = r y/2—(3/N. 
The density picture correctly describes the BA roots wherever <7 C (A) is of order 1. 
From eq.(4.18) we then find as validity interval 

\\\<0{^JjN) (4.19) 



which shrinks to zero when iV — > 00. Roots \Xj\ > 0(y//3/N) have a spacing of 
order larger than 0(1/N) and cannot be described by densities. In particular, the 
roots with largest magnitudes have finite iV — > 00 limits spread by 0(1) intervals 
(we checked this fact numerically too). Therefore, contrary to the usual situation 
[15], we must go beyond the density description to obtain a bulk quantity like the 
free energy per site. 

We now consider the limit iV — > 00; iV explicitly enters eqs.(4.16) and (2.25) 
only through the functions iPn(X) and z c (X). Thus, inserting eqs. (4.16) and (4.17) 
in eq.(4.14) yields 

f(/3, h) = E xxz + p-'HP, h) (4.20) 

where Exxz is the ground-state energy at zero external field h of (4.1), namely 

+00 



Exxz = 2J 



. sinh(7r — 7) k 

cos 7 — sin 7 / dk 



(4.2i; 



sinh7r/ccosli7A; 
while (the spin reflection symmetry allows to take the half sum of the two cases 
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h > and h < 0, to ensure explicit reality of the final expression) 



-hOO 

L(8h)-1m /d\ l0g[1 + e ' £+(A)][1 + e ' £ ' (A)] (4 221 

L(/3, h) - Im y rfA 27sinh[7r(A + ^ )/7] ( 4 - 22 ) 



The functions e±(A) satisfy the limit N — ► oo of eq.(2.29), which stays formally 
unchanged 

e ± = g ± -G* log(l + e~ €± ) + G 2v * log(l + e" e >) (4.23) 

but where now 

2nij3 ir3h 

7 sinn[7rA/7j 7r — 7 

We may also write e±(A) = — «Z(A + 277; =f/i) (compare with eq.(2.28)), with the 
analytic function Z{\) simply related to the original counting function by 



Z{\)= lim [Z jV (A)-7rATsign(A)] (A real) (4.25) 



The problem of calculating the free energy of the (infinite) XXZ spin chain has thus 
been reduced to the problem of solving the two coupled nonlinear integral equations 
(4.23) and performing the integral (4.22). Notice that for vanishing as well as for 
purely imaginary external field, eq.(2.28) entails e + = e_, so that eq.(4.23) reduces 
to a single equation for e(A) = — iZ{\ + irj). 

Notice that for A — >■ ±00, eqs.(2.11) and (4.25) yield 

Z(±oo) = -2u 

One could take the limit iV — >■ 00 directly in the definition itself of the counting 
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function, eq.(2.11), which now reads 

N 
Z N (\) = S N {\) - J2 0(A - A fc , 7) + 2i/3h (4.26) 

fe=i 

with the source term 

S N ( A) = N U( X + i-f/2 - ifi/N, 7/2) + <f>( A - * 7 /2 + i/3/JV, 7 /2)l (4.27) 

Since the roots Xj all have finite A" — > 00 limits and 

lim [SW(A) - tt Assign (A) - N<f>(\, 7)] = 2/3g(A) (4.28) 



with 



sinh 2A , , 

g(A) = — : — coth A 4.29 

yv ; cosh 2A- cos 2 7 v ; 



we obtain for the limit (4.25) of the counting function 






Z{X) = 2(5q{X) - ^ [0(A - 0,7) - 0(A, 7)] + 2i(3h (4.30) 

jell 

where the real numbers £j, which satisfy Z(£j) = (2j — 1)tt, j G 7L, are the 
iV — > 00 limit of the original BA roots Ai, A2, • • • , Ajv- Eq.(4.30) shows that Z(X) 
has periodicity in and has, as unique singularity in the strip |Im A| < 7/2, a simple 
pole at the origin with residue —2/3. This fact will play an important role in the 
high temperature expansion of sec. 6.1. 

The new BA equations Z(£j) = (2j — l)ir embody all the information about the 
XXZ thermodynamics. They form an infinite set of algebraic equations equivalent 
to our non-linear integral equations (4.23). We want to stress that these roots are 
discrete although we have already set A^ = 00. The root spacing is actually of 
order T and in the T — > limit we recover the continuous distribution of roots 
associated with the antiferromagnetic ground-state [15]. The roots £j have an 
accumulation point at £ = 0. That is, £j = — (3/[k(j — 1/2)] for j — ► 00. For low 
temperatures, the largest root, £1 , is of order (7/vr) ln/3. 

26 



We can verify that eq.(4.30) is indeed equivalent to eq.(4.23) by means of the 
usual contour integrations. Taking into account that Z(X) has a simple pole at the 
origin and using the identity 

q(X) = -i<t>'{\ - z 7 /2, 7 /2) + \i<t>'{\ 7) (4.31) 

one obtains from eq.(4.30) , for 7/2 > ImA > and h > 0, 

Z(X) = 2i(3 <t>'{\ - z 7 /2, 7 /2) + 2i(3h - I ^L. <t>'{\ - /x, 7) log[l + e~ iZ ^\ (4.32) 

J 2m 
r 

where the contour V encircles all the roots £j as well as the origin. By deforming T 
into straight lines as in sec. 2, applying the convolution (1 + K)^ 1 and observing 
that 



+00 

d\x (1 + K)~\\ - ix) ${\x - i-i/2 + iO, 7/2) 



7sinh[7r(A + i0)/7] 

( 4 - 33 ) 



dfx(l + K)- L (X-fj) 
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2(tt- 7 ) 



one recovers the N — »• 00 limit of equation (2.25) (with u = —i(3h and 6 = 
7/2 — (3/N), from which eq.(4.23) has been derived. 

As we shall see in sec. 6.1, the alternative form (4.32) of the integral equation 
for Z(\) is most suitable to study the high temperature regime, since it involves a 
complete contour integral. It is then useful to obtain an analogous formula for the 
free energy. From eq.(4.22) we obtain (for h > 0) 



+00 



1 r Z(X + ^0) 1 r^ 



L(/3, h) = - dX , ,. v - ' — ^— n - — (bdX L , ,. — J (4.34) 



1 + e ^ z ( A ) 



2 7 7sinh[7r(A + i0)/7] 2i J 7sinh[7rA/7] 

-00 ' r 
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Substituting eq.(4.32) for the term linear in Z{X) and using (4.33) , leads to 



^ ~ I « * l \-JZ&^ +m +f£<^MW~«»t 



By Fourier transformation, the first integral can be rewritten as 



(4.35) 



+00 



d\ , = /3 (2 J cos 7 - E X xz + /i) (4.36) 

7Sinn[7r(A + 2OJ/7J 



so that we finally obtain using eq.(4.20), (4.35) and (4.36) 



f(P,h) = h + 2Jcosj + 



d\ lo § 



1 + e- iZ W 



sin 7 

/3 7* 27ri sinh A sinh(A — 17) 

r 



(4.37) 



It is clear that, to go from eq.(4.22) to eq.(4.35), we have just performed in reverse 
order the same operations that led to (4.22) from the definitions (4.14) and (4.15). 
However, in the middle, we have taken the limit iV — ► 00 , which is difficult to 
take directly on eqs.(4.14) and (4.15) . The final result, eq.(4.37), can now be 
transformed, via contour integration, into an infinite sum over the roots £j with 
the result 



f(P) = h + 2Jcos 7 + -^ logj]! 1 + 2 cos 7 ) 

P k 



2sin7sin(sfc — 27) 
cosh(rfc) — cos(sfc — 37) 



where r^ = 2 Re (£&) and Sk = 2 Im (£&). 



28 



5. Ground state scaling function of the 
massive Thirring (sine-Gordon) model 

Let us now reinterpret the diagonal-to-diagonal lattice associated to the trans- 
fer matrix i2i\r(A, 0,ou) as the unit-time evolution operator for Minkowski space- 
time discretized in light-cone coordinates. That is, the axis correspond to x ± t, 
(x and t being the usual space and time variables). For this purpose we choose 
= — 2iO, with © real. Then, the local R— matrices (2.2) are rendered unitary 
upon multiplication by o(— 2i@) -1 , 

In this light-cone approach, we start from the discretized Minkowski 2D space- 
time formed by this regular diagonal lattice of right-oriented and left-oriented 
straight lines. These represent true world-lines of "bare" objects (pseudo-particles) 
which are thus naturally divided in left- and right-movers. The right-movers have 
all the same positive rapidity 0, while the left-movers have rapidity — 0. One can 
regard as a cut-off rapidity, which will be appropriately taken to infinity in the 
continuum limit. Furthermore, we shall denote by V the Hilbert space of states of 
a pseudo-particle (we restrict here to the case in which V is the same for both left— 
and right-movers and has finite dimension n — 2, although more general situations 
can be considered [20,15] ). 

The dynamics of the model is fixed by the microscopic transition amplitudes 
attached to each intersection of a left- and a right-mover, that is to each vertex 
of the lattice. This amplitudes can be collected into linear operators Rij, the local 
R— matrices, acting non-trivially only on the space V{ ® Vj of ith and j'th pseudo- 
particles. R{j thus represent the relativistic scatterings of left-movers on right- 
movers and depend on the rapidity difference — (—0) = 20, which is constant 
throughout the lattice. Moreover, by space-time translation invariance any other 
parametric dependence of Rij must be the same for all vertices. We see therefore 
that attached to each vertex there is a matrix i?(20)^, where a, b, c, d are labels 
for the states of the pseudo-particles on the four links stemming out of the vertex, 
and take therefore n distinct values. This is the general framework of a vertex 
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model. These local R— matrices (2.2) were rendered unitary upon multiplication 
by a(-2i0)- 1 . 

The difference with the standard statistical interpretation is that the Boltz- 
mann 

weights are in general complex, since we should require the unitarity of the ma- 
trix R. In any case, the integrability of the model is guaranteed whenever i?(A)^g 
satisfy the Yang-Baxter equations. 

For periodic boundary conditions, the one-step light-cone evolution operators 
Ul(Q) and Ur(Q), which act on the "bare" space of states Hn = (®V) 2N , (JV is 
the number of sites on a row of the lattice, that is the number of diagonal lines), 
are built from the local R— matrices i?y as [13] . 

Ur(&) = W(G)V , U L (G) = W^V- 1 , W(G) = R 12 Rs4 ■ ■ ■ R2N-1 2N 

(5.1) 
where V is the one-step space translation to the right. Ur ( Ul ) evolves states 
by one step in right (left) light-cone direction. Ur and Ul commute and their 
product U = Ur Ul is the unit time evolution operator. 

If 5 stands for the lattice spacing, the lattice hamiltonian H and total momen- 
tum P are naturally defined through 

U = e~ iSH , Ur Ul 1 = e iSP (5.2) 

These operators directly express in terms of the alternating transfer matrix £2jv(A, — ^©, 0)| 
as follows 

U R (&) =a(-2i&)- N t 2N (-i&,-i&,0) , U L (&) = a(2i&) N t 2N (+i®, -i®, 0)" 1 

Therefore, the unit time evolution operator results 

t(-iO, -iB, 0) t(iO, -iG, 0)" 1 (5.3) 



U = e~ im 



a{2iQ) ' V 



_a(-2iQ)_ 

Notice that H and P are (non-local) lattice operators. For simplicity we set the 
twist 10 = in this section. Whenever the R— matrix, that is the building block of 
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the whole construction, satisfies the Yang-Baxter equations, this UV-regularized 
QFT is integrable, and will therefore be so also in the continuum limit. The 
integrable QFT corresponding to the six-vertex R— matrix given above was iden- 
tified in ref. [13] by Jordan-Wigner transforming the local Pauli matrices into 
discretized fermionic fields (on the infinite lattice), and then showing that these 
fields satisfy a lattice version of the field equation characteristic of the massive 
Thirring model. The bare continuum fields are recovered when S — > and simul- 
taneously — ► oo with the bare mass scale mo ~ <5 _1 exp(— 20) held fixed. The 
renormalized , physical continuum limit follows instead by keeping fixed the true 
mass scale m ~ 5 _1 exp(— irQ/^f) (0 plays essentially the role of cutoff in rapidity 
space). All the on-shell calculations can be performed exactly within the Algebraic 
Bethe Ansatz, leading to the complete spectrum and exact S— matrix of the mT 
model. Let us also recall that the latter is equivalent a la Coleman-Mandelstam 
to the sine-Gordon model. We shall now show what modifications to the general 
formulae derived in sect. 4 are necessary in order to discuss the thermodynamics of 
the sG-mT model. 

By the euclidean symmetry of the corresponding functional integral, we know 
that the free energy density of a 2D QFT at temperature T (on an infinite line) is 
identical to T times the ground state energy of the same QFT on a ring of circum- 
ference /3 = 1/T. Of course, such a quantity is UV divergent and must be properly 
subtracted. Requiring that the free energy f((3) vanishes at zero temperature is 
equivalent to consider only the Casimir energy, that is the difference between the 
ground state energy on the circle and that of the infinite line. In the context of 
Perturbed Conformal Field Theory, this Casimir energy is known as ground-state 
scaling function E(/3), so that we have f(/3) = E(f3). When T — ► oo we expect 
that UV fixed point characterizing the QFT manifests itself in some universal fash- 
ion. In detail, we expect that E(/3) tends to the Casimir energy of the Conformal 
Field Theory corresponding to that UV fixed point. From this point of view, the 
situation of a massive QFT such as the sG-mT model is opposite to that of a 
gapless spin chain such as the XXZ model at < 7 < n. In the former the mass 
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is generated by some IR relevant perturbation of the CFT and becomes irrelevant 
in the UV limit T — ► oo. In the latter the lattice cutoff corresponds to the pres- 
ence of IR irrelevant perturbations of the CFT those effects become negligible in 
the T — > limit. Thus the basic properties (central charge and scaling dimen- 
sions) of the same CFT (the free massless boson) will appear at low temperatures 
in the gapless XXZ spin chain and at high temperatures in the massive sG-mT 
model. This dichotomy will be quite clear in the unified description that we are 
here proposing, since, apart from a different source term, the fundamental integral 
equation is essentially the same in both cases. 

As stated above, the diagonal-to-diagonal transfer matrix T/v(0), eq.(4.2), 
once unitarized into U according to eq.(5.1), describes the real time evolution of a 
discretized mT model. The continuum relativistic QFT, on the infinite Minkowski 
plane, is obtained taking first the IR limit iV — *■ oo at fixed lattice spacing S and 
then the continuum limit S — » near the critical point = oo (with the physical 
mass scale m ~ 5 _1 exp(— irQ/'j) held fixed). On the other hand, by taking the 
continuum limit at fixed NS = /3, we get instead the same QFT on a ring of length 
/3. In particular, if we consider the ground state on the lattice, this limit will 
yield the (UV divergent) bulk ground state energy of order (3 plus the ground state 
scaling function E(/3). 

Using now eq.(3.4) and (5.3), the lowest value of the energy is then attained 
in the a.g.s. configuration and take the form 



N 

En = J2 M X 3 - > 7/2) - 0(A, + 6, 7 /2) - 2vr] (5.4) 



where the BAE roots Xj are real, since u = 0. In terms of the counting function 
Zjv(A), the expression for E^ is identical in structure to eq.(4.16), with only E c 
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and ^jv(A) changing to the forms proper of the mT model: 



N 2 

Er= 

(3 

N 



-2tt + / dX 



( KA + 29,7/2) 
7cosh7rA/7 



MX) = —q [sech f (A - O) - sech ^(A + 0) 



(5.5) 



In the continuum limit N — > oo the rapidity cutoff diverges like 



~ _ log — - 

7i rap 



(5.6) 



where m is the mass of the mT fermion (or sG soliton). Thus 



E c = jTti f3 cot |- + UV divergent terms 



(5.7) 



where the first finite term is the scaling bulk energy (the same result was obtained 
in ref.[19] by completely different means). By 'UV divergent terms' we mean terms 
that scale as the bare mass (e ) for large 0. On the other hand 



E(0) = lim [E N - E c 

+00 



- — Im / dX sinh [tt(X + if))/"/] log 

1 J 

— oo 



l + e 



iZ(\+ir}) 



(5.1 



where Z(X), the A^ — > oo limit of the counting function Zjy(A), satisfies 



+oo 



Z(X) = m(3 sinh(7rA/7) + 2 Im dpi G(X — /J, — irj) log 



l + e 



iZ(fi+iri) 



(5.9) 



Let us recall that this result applies to the complete range < 7 < n, provided 
we take < 77 < 2 m i n (7? ^ — l)- in the repulsive regime 7 < n/2 we can cast 
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eqs. (5.8) and (5.9) in a very nice form, reminiscent of the standard TBA. Let us 
choose i] = 7/2 and set 

6,(5) = -»Z( Je + < 7 /2), e f = TJ 

L f = log(l + e e , Lj = log(l + e e /) = I7 (5.10) 

G o (0) = 2G(J0) , Gi(0) = G o (0 + m - iO) 

Then the ground state scaling function reads 

+00 
m f r 

E(P) = / dO cosh^ L f (6) + Lj{6) (5.11) 

— CX3 

while the nonlinear integral equation becomes 

€f = g-G *L f + G 1 *Lj (5.12) 

where now g(8) = m(3 cosh 6 is the source term of the standard TBA, and is directly 
related to the level density of free relativistic particles. Moreover, from eqs. (2.24) 
and (5.10), 



47r sinh (f^ - i) A; cosh f A; 



By construction, the equation for ei is just the complex conjugate of (5.12). Let us 
observe that Gq(9) is related in the standard way to the fermion-fermion scattering 
amplitude Sff(9) 

G ^ = h^ oiS " w (5 - 13) 

while, up to a sign, G\{9) coincides with the crossed Go, that is Gq{itx — 9). In 
other words, we can regard ej and ej as complex pseudoenergies for the fermions 
and antifermions, respectively, with eq.(5.12) as TBA equations involving only 
physical particles. This situation should be compared with the standard TBA 
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for the mT (or sG) model [23], which for nonrational values of '-f/ir involves an 
infinite number of nonlinear integral equations for the infinitely many different 
types of "magnons" describing the energy degeneracies of the physical fermions and 
antifermions. Therefore, our equation effectively provides a resummation 
of the magnon degrees of freedom. 

The results (5.11) and (5.12) may be extended to more general situations in a 
rather straightforward way. For instance, we could consider a nonzero twist uj in 
eqs.(5.3), or introduce a coupling to the conserved U(l) charge Q in the real time 
evolution defined by U. Let us consider first this latter case. The change of the 
Hamiltonian H — ► H — hQ, (h > 0), has a simple consequence: for h > m the 
ground state contains a definite number of fermions (which have charge Q — 1) 
with the lowest possible kinetic energy. On the light-cone lattice this means that 
the solution of the BAE must contain holes around A = 0, that is roots £j of the 
counting function Zjy(A), Zjv(£j) = which do not appear in the sum over BAE 
roots in the definition of Zjy(A). 

As in the XXZ chain, we can find for the sG-mT model an infinite set of alge- 
braic Bethe Ansatz type equations which are equivalent to the non-linear integral 
equation (5.9). In order to do that, we define 



v(6) = Z{~f0/n) - mpsmh(9) 



(5.14) 



This function v{9) vanishes at 9 = ±oo. We can then write, 



Im log 



1 + e 



iZijO/ir+iO) 



-v(0)+«D(0) 



(5.15) 



where D{9) is the (subtracted) discontinuity of log 



1 + e iZ(ye/ir+iO) 



. That is, 



I 1 4. e iZ(~/8/-K+iO) a 

Diff) = _ log tJLl — — - + i!2i sinh(^) 



(5.16) 
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Then, 



D '(0) = V± cos h(0) - J2 W - k ) 

k 



(5.17) 



where the Q k fulfil the equations: 



Z("f6 k /n) = (2k + l)n, keZ 



(5.18) 



Since Z(— t 
origin: 



-Z{6), the 6 k are symmetrically distributed with respect to the 



0fc = — 0-k-l 



(5.19) 



Using the asymptotic behaviour Z(j9/tt) = m(5 sinh(0) , we find for large roots 
Qk , 



Ok - ± log 



2tt 
m(3 



|2fc + l| 



(5.20) 



We see from eq.(5.20) that the term in cosh0 in eq.(5.17) simply ensures the 
finiteness of D'(6). As eq.(5.20) shows, the roots accumulate at infinity in the mT- 
sG model, whereas they accumulate at the origin for the XXZ chain at temperature 
T (sec.4). 

Eq.(5.9) for rj — »■ + can be now written as 



+ 0O 



1/(0) = dfi G (6 - fj) [u(fj) + 2ttD(ii)] 

—00 

Upon Fourier transform and pulling i/{9) to the l.h.s., this yields 

+00 

K0) = (z) 2 [dn<t>'((l) 2 [e-n], 1 )D(n) 



(5.21; 



7T 



7T 



(5.22) 



or 



7 N 2 



Z( 7 0/7t) = m(3 sinh 6 + dfi </>((-) 2 [0 - ju] , 7) £>'(/x) 



7T 



(5.23) 



This is the analog of eq.(4.32) . Inserting now eq.(5.17) into eq.(5.23) we find for 
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7 > 7r/2 

Z( 7 0/7r) = -J>((^) 2 [0-0 fc ], 7 ) (5.24) 



where the sum is to be understood in principal part. This shows explicitly that 

2 

Z(X) has ^- as period. Moreover, setting 6 = Q\ in eq.(5.24) and using eq.(5.18) , 



T 2 

7 

we find 



(2/ + ]> = - J2 0((-) 2 [9i ~ k ] , 7 ) (5.25) 

k 
always for 7 > tt/2. 

6. The Recursive Regime 

Depending on the temperature regime chosen, the integral equations (4.32) 
and (5.12) can be simply solved by iterating the inhomogeneity. This happens for 
small j3 in the XXZ chain [eq.(4.32)] and for large (3 in field theory [eq.(5.12) for 
the niTm ]. We call 'recursive regimes' those where such expansions apply. 

6.1. High Temperatures in the XXZ chain 

We study here the free energy f{f3) of the XXZ chain for high temperatures. 
When /3 is small it is convenient to use the form (4.32) plus the uniform expansion 

00 

Z(\) = J2P k h(X) (6.1) 

fe=i 

Since the residue of Z(X) is linear in (3 only , 61(A) has a pole at the origin, we 
have 

6 1 (A) A ~°-^ + 0(A) (6.2) 

while all the 6fc(A) for k > 2 are analytic there. 
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From eq.(6.1) we find up to fourth order in j3 



log 



1 + «r**W 



ln2-| 61(A) 



2 

2 



MX) + i &i(A) i 



2 



^s(A) + - 61(A) 62(A) 



164(A) + i 61(A) 63(A) + ^6 2 (A) 2 + 1&!(A) 4 



0(^ 5 ) 
(6.3) 
We insert this expansion into eq.(4.32) and notice that only the poles at A = in 
eq.(6.3) contribute to the contour integral over T . Since only 61(A) is singular, the 
procedure is perfectly recursive and we find 



61(A) = 2q(X) + 2ih 
62(A) = -l 



</>"(A,7) -ih4\\l) 



63(A) = ih(f)'(X, 7) cot 7 



64(A) = — ih 

ih 



h 2 



., . 0'(A,7) + 

6 2sm 7 3 / 3 V 3 



* +VU'(A, 7 ) 



sin 7 



^0 W (A,7)-— W/ (A,7) 

(6.4) 
where 6 = /i/(J sin 7). Next we insert the expansion (6.3), with the explicit results 
listed in (6.4), into the formula (4.37) of the free energy. Again all we need to do 
is to repeatedly apply the residue theorem, with the final result 

/(/?) = - p- 1 In 2 + J cos 7 - (3 [ J 2 (l + ±cos 2 7) + \h 2 ] + (3 2 cos 7 ( J 3 + 6 2 ) 

- /9 3 [(i + I cos 2 7 - A cos 4 7)^ - J 2 6 2 - ±h A ] + 0(/5 4 ) 

(6.5) 

which indeed agrees with the high T expansion, as can be derived direclty from 

the definitions (4.4) and (4.1). 
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6.2. Low Temperatures in field theory 



For low temperatures (large mf3) the non-linear integral equation (5.12) valid 
for 7 < 7r/2, 



e/(A) = m/3coshA + 

+oo 

dji \ G\{\ — fx) log 



1 + e 



-e/(^) 



Go(A-//)log 



1 + e -e 'W 



(6.6) 



can be easily iterated using the inhomogeneity mj3 cosh A as zeroth order approxi- 
mation. That is, 

c/(A) = /o(A) + /i(A) + / 2 (A) + ... 

/o(A) = m/3cosh A 



+CX) 



/i(A) = I dti[ Gi(A - fi) - G (A - //) ] log 
dfi 



1 + e 



-m/3cosh^i 



(6.7) 



-oo 

+00 



/2(A) 



=/o(m) -f l 



Go (A - /i) fi(fj.) - Gi(A - /i) /i(/x) 



The energy i?(/9) is expanded in analogous way as a sum: 



£03)=£o(/3) + £i(/3) + £ 2 (/3) + .. 



(6i 



where 



+oo 



111 



E>o(P) = / dfi cosh/i log 



1 + e' 



-m/3 cosh/i 



— oo 
-00 



(6.9) 



E ^ = ^ I dficosh^e-^) /l(/i) + /l ( /x ) 



Let us now analyze more explicitly Eq(/3) and E\(f3). From eq.(6.9) we find 



2m 



w / i\n— 1 

W) - -— J2 — Ki(nm/3) (6.10) 

n—\ 

This is the typical behaviour of a massive QFT. Notice that eq.(6.10) for 7 = 7r/2 
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gives the exact energy E(/3). (In that free case Gq(A) = 
For large m/3, eq.(6.10) yields 



E Q {(5) P =°° -—K 1 (mP) + 0(e~ 2m P) * 



1m 



-m(3 



1 + Oi 



a' 



-0( e -2m/3) 
(6.11) 



It is easy to see that E n {(3) is of order e \ n + l ) m P for large m(3. 
We get from eqs.(6.7) and (6.9) 



+CX) 



+oo 



/;/ 



E 1 {(3)= --^ / dX coshAe- m/3coshA / d/x log 

— 00 — OO 

[ 2 G (A - jk) - Gi(A - ix) - Gi(A - /x) } 



1 + e 



-m/3 cosh /i 



(6.12) 



The terms inside braces [. . .] in eq.(6.12) admit for 7 < tt/2 an integral represen- 
tation that follows from eq.(2.24) and eq.(5.10) 

+00 
2 G (A) - Gi(A) - Gi(A) = - | ^ e^ A [1 - ^)] (6.13) 



where 



sinh(4-^)|/ + 3sinh(^-2)y 
2 coshy sinh(- — l)y 



(6.14) 



For large m/3 we can approximate the log in eq.(6.12) by the first term in its power 
expansion. This yields, 



+00 +00 +00 

E 1 ((3)=-— ! d\ coshA e- m/3coshA / d/x e - m(5cosh ^ ( \ 



+ 0(e- 3m(5 ) 



e ix( - x ~^ [1 - of— 11 
2tt [ yv 2 ;J 



(6.15) 



This integral can be approximated by the saddle point method for large mf3 with 
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the following result: 



El (P) P =°° e- 2m ? 



m K y _ , 1 



27T/3 (3 



02' 



+ 0(e~ 3m/3 ) (6.16) 



where 



K^ 



dy i v 

— 9 2/ 

7T Z 



+CX) 



rf y sinh(4-^ + 3sinh(^-2)y 
7r 2 2 coshy sinh(^ — l)y 



(6.17) 



This function of 7 can be expressed as an infinite sum of ip functions. For rational 
values of 7/71 it admits simpler expressions. For example, 



tfo = -i(41n2-l) 

7T Z 



K, 



tt/3 



7T 



In addition, we see from eq.(6.17) that K^ < for 7 < n/2 and that lim 7 ^ 7I -/ 2 - K^ =1 
—00. This divergence indicates a change of regime since E\(f3) is in fact identically 
zero at 7 = ir/2. 

In summary, E n (/3) for large m/3 is at leading order 0(e~^ n+1 ' m ^) and contains 
contributions with arbitrary higher powers of e~ m ^. 

From Eq(/3) + E\((3) we can write E(/3) up to contributions 0(e~ 3m ^). We 
find from eqs. (6.16) and (6.11) , 



7T y Tip 



1 + V2 



^ + °$ 



-0(e 



-3m/3\ 



The full term of order e ^ in E\ {(3) possess the following integral representation: 

+00 
El ((3) P =°°- [ dx K ix {m{3) [K 1+ix (mf3) + K^^m/3)] 

71" ./ 



sinh (irx/2) sinh (nx[-^- — 1] 
cosh(7ra/2) sinh (vrx[^ — 1] 



+ 0(e 



-3m/3\ 



where we used eqs. (6. 12), (6.13) and (6.14) . 
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7. The conformal regime 

The conformal regime is the opposite to the recursive regime. That is, large f$ 
for magnetic chains and small /3 in field theory. 

In the conformal regime the calculations are more subtle but still straightfor- 
ward as shown below. 

The starting point is again the integral equation (2.25). For simplicity we shall 
consider here only real twists u, that is purely imaginary external fields h for the 
XXZ chain. Then, Z(X) = Z(X) and we can write for all cases 

+00 
-ilogF(fi) = v?(/i) + 2Im d/j,' G(fi- fi' -ir])\og[l + F(/j,' + ir])] (7.1) 

—00 

where F(/j) = e iZ W anc [ 



sinh(7r^/7) 
cos(7r6>/[27]) 



7T — 7 



finite size vertex model 



2N arctan 

<p(fj) = < . y + H^ XXZ thermodynamics ( 7 - 2 ) 

m(3 sinh(7r/i/7) sine-Gordon field theory 

The solution Z(X) must be then inserted in eq.(3.9). [ Recall that e(/i) = —iZ(fi)] 



7.1. Leading finite size corrections for the six- vertex model 

Let us derive here the dominant finite size corrections for the six-vertex model 
free energy [see sec. 3]. 

We see from eq. (7.2) that </?(//) ~ iV for large iV as long as /j, < (7/V) log[iVcos(^-)].l 
Hence, Z(/j,) ~ N and we find that this bulk contribution to L^ vanishes to the 
order iV° and iV -1 . 
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The contributions to Ln of order N (and smaller) come from values of |/x| 
larger than (7/vr) log[iVcos(— )], where the previous estimate Z(fi) ~ N does not 
hold anymore. In order to compute the contributions from large positive values of 
fi one introduces the new function 



F(x) = e'<^ , x = fi-- \og[AN cos(— )] (7.3) 

7T 7 



Then eqs. (7.1) and (3.9) reduce, in the N — *■ 00 limit, to 



+00 

TTX 2,TTC0 f 

—i\ogF(x) = —e y Im / dy G(x — y — irj) log [1 + F(y + irj)] (7.4) 

7T - 7 J 

— 00 



L 2A r = iw - — -—3- / rfx e 'Im e 7 log[l + F(x + irj)} \ (7.5) 

A^7 cos(^^) J < > 

' — CXD 

For simplicity, we shall choose rj = + and use that G(x+iO) = G(x) and G(±oo) = 
for real x. We shall analogously compute below the contributions from large 
negatives values of //. 

Fortunately, it is not necessary to solve the integral equation (7.4) in order to 
calculate the integral (7.5). In fact, we can use the following lemma 

Lemma. Assume that F(x) satisfies the nonlinear integral equation 

-i log F(x) = <p(x) + 2 / dy G(x - y) Im log[l + F(y + iO)] (7.6) 

Xl 

where tp(x) is real for real x and x\, X2 are real numbers. Then the following 
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relation holds 



Im /*•(,) k«[l + F(* + <0)] = 

XI 

\ Im [<p(x 2 ) log(l + F 2 ) - ¥»(xi) log(l + Fi)] + \ Re [£(F) - £(F 2 )] 

+ / dj/ [ C(x 2 - y) log(l + F 2 ) - G(xi - y) log(l + Fi) ] Im log[l + F(y + iO)} 



Xl 



(7.7) 



where Fi j2 = F(xi ]2 ) and £(t) is a dilogarithm function 



e(t) = du 





log(l + -u) log-u 
■u 1 + u 



(7- 



To prove the lemma we consider the relation 



X-2 



(( F 2 )-({ /-j ) = / <M log[l + F(x + iO)]—\ogF(x) - log F(x) — log[l + F(x + z( 



.-f'l 



and then use eq.(7.6) and its derivative to substitute logF(x) and d log F(x)/dx. 
This yields 



J'2 



£(F 2 ) - £(Fx) 



; / ( ] A . <; J^5> i g[l + F(x + iO)] - y»(x) — log[l + F(s + iO)l 

dx dx 



•fi 



X'2 

+ log[l + F(s + iO)] / dy G'(x - y) Im log[l + F(y + i( 



.Tj 



— log[l + F(x + iO)] / dy G(x - y) Im log[l + F(y + iO)] 



Xl 



The double integral here cancels upon partial integration and using that G'(x—y) = 
—G'(y — x). Finally, taking real part yields eq.(7.7) (related identities were used 
in ref. [18]). 
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We can now use the lemma to compute the integral (7.5) setting 



, , _ZI£ 7100 

ip(x) = — e i and x\ = — oo , X2 = +oo. 

tt — 7 



We have F(xi) = 0, F(x 2 ) = e~ ZluJ and 



+CX) 



71 



7 



— dx e i log[l + F(x + i( 



- -M 



v. 2 ^) + £( e +^ 



UJ 



2 1 



'7T 



(7.9) 



In this expression we can apply the formula 



e(z)+e(i/z) 



7Y^ 

y 



(7.10) 



The contribution from large negatives values of \x is computed analogously intro- 
ducing the variable 



7 



' i ' i 

x = ji H — lo 



TX 



4Ncos(^) 



We find from eqs. (7.5) ,(7.9) and (7.10) collecting both contributions : 



71 ( 7\\ 

L2N = 6iV tan {- 



Qu z 



2(l- 7 /vr) 



7T 



(7.H) 



Hence the central charge turns to be unit. The factor tan ( — ) is a geometric effect 
due to the kind of diagonal-to-diagonal lattice we are using. (In the row-to-row 
framework a factor sin ( — j arises [15]). 

Let us prove that this tan factor is just the speed of sound for low-lying exci- 
tations. The eigenvalue of logi2iv(A, 0,uo) for a hole excitation normalized to the 
a.g.s. writes [15] 

(7.12) 



— 1i arctan 



e ^(0+iA) 



where stands for the position of the hole. Since we are here on an euclidean lattice, 
we can write log£2iv(A, 0,oo) = h + ip to define the time and space generators h 



45 



and p, respectively. For large positive or negative 0, the eigenvalues of h and p 
result from eq.(7.12) 

e = ±2 sin — e i , p — ±2 cos — e i 
\1 J V 7 / 

That is, the speed of sound turns out to be precisely tan ( — ) • 

This result can be interpreted as a proof of the presence of conformal invariance. 

7.2. LOW TEMPERATURES IN THE XXZ CHAIN 

We shall consider now the low temperature regime. When (3 3> 1 eq. (4.23) 
indicates that Z(X) ~ /3, so that log[l + e ' '] ~ iZ(X), at least as long as 
|A| < (7/vr) log/3. At dominant /3 ^> 1 order eq.(4.23) linearizes in the same 
way as it does for small (3. Inserting then Z(X) ~ 2f3q(X) in eq.(4.22), yields 
lim^-KXj L(/3) = 0. Therefore eq.(4.20) tells us that 

fm^Exxz^ + Otp- 2 ) (7.13) 

The contributions of order (3~ 2 (and smaller) come from values of A larger than 
(7/71) log/9, where the previous assumption log a ~ O{0) does not hold anymore. 
It is then convenient to introduce the new function 

F{x) = e iZ ^ , x = A - 1 log ^ (7. 14) 

71 7 

Then eqs. (4.22) and (4.23) reduce, in the (3 — > 00 limit, to 



r/3 
and 



+00 
L(/3) = i /"rfa;e^ x/7 Imlog[l + F(a; + iO)] 

7T/9 i 



(7.15) 



+00 



-i log F(x) = -e~ ra/7 + 2 f dy G(x - y)Im \og[l + F(y + i0)] (7. 16) 

— 00 

The integral in eq.(7.15) may now be exactly calculated by using the lemma 

4G 



from section 7.1 with (p = — exp(—7rx/'~f) and x\ = — oo, x 2 = +00. We have 
F(xi) = 0, F(x 2 ) = 1 and 



+CX) 



21m [ dxe-™/t log[l + F(x + iO)] = -^(l) = -^ 

J 7T 6 



Then the free energy for low temperature reads 

f{P) = Exxzd) - J2f-— P~ 2 + °(P~ 2 ) ( 7 - 17 ) 

in perfect agreement with refs. [2,16,17]. 

As we have just shown, both the high and the low temperature leading be- 
haviors of the free energy can be derived without much effort from our non-linear 
integral equation (2.27). The higher order corrections for high T can be obtained in 
a very systematic way. The situation for the o(/3 -2 ) terms in the low T expansion 
(7.17) is more involved. 

However, the NLIE (4.23) admits a Riemann-Hilbert formulation analogous to 
the one presented in sec. 7.4 for the sG-mT model. One can establish from such 
Riemann-Hilbert formulation the analytic structure of the low T expansion. 

7.3. High temperatures in the mT-sG field theory 
Let us recall the integral equation (5.12) 

e = rch-Go*L + Gi*l (7.18) 

where e = ej, r = m/3 and ch stands for the hyperbolic cosine, ch($) = cosh#. 

For small r (high temperatures), the regions with positive and negative where 
rcosh# ~ 1 are very far apart. It is then useful to separately study eq.(5.12) with 
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inhomogeneities 



r _ ±6 _ ±e-log2/r 

2 



That is, we define the kink function e k as the solution of the same integral equation 
with r ch replaced by an exponential 

e k = e -Go*L k + G 1 *L k (7.19) 

where L k = log(l + e _£fc ). (Analogous kink functions were first considered in [4] 
and have now become standard in TBA calculations). 

Then, we set 

7j(0) = e{9) - e k (9 - log2/r) - e k (-9 - log2/r) 

(7.20) 
£{9) = L{9) -L k (9- log2/r) - L k (-9 - log2/r) + log 2 

From eqs. (7.18) and (7.19), using the translational invariance of eq.(7.19) 77 and 
£ turn to be related by an homogeneous equation 

r] = -G *£ + G 1 *£ (7.21) 

We then find for the energy 

+00 
777 / — 

E(P) = / d9 cosh# \L{9) + L(9)] 

2-k J 

(7.22) 

+00 v ' 

= "^ (/* + /*) " ^ / & cosh 9 W) + 1 ( )} 
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where the kink contribution reads 

+00 



f k = J d9e e [L k (9 - log2/r) + L k {-9 - log2/r) - log 2] 

— OC 

+OO +OO 

= 2 - Jd9e e L k (9) + J d9 e- e ^[L k (6 -log2/r) -log 2] (7.23) 



-00 
+00 



+00 



lf d0e » Lm + lf d e e -» 8 J± 



The real parts of these two integrals can be calculated without explicit knowledge 
of the kink function e(9). To the term proportional to r _1 we can apply the Lemma 
with the substitutions f{6) = e and x\ = —00, X2 = +00. This gives 



^00 

Jd9e e [L k (9) + L k (9)} =^(l) = y 



(7.24) 



To compute the term linear in r in eq.(7.23) we can adapt to our situation an 
argument in ref. [4] : to converge for 9 — ► —00, the second integral requires that 
dL k /d9 decays faster than e , hence also e k (9) must tend to faster than e as 
9 — ► —00. We then differentiate eq.(7.19) with respect to 9, let 9 — ► —00 and use 
the asymptotic expansion for large |A| 

G (A) ~ 01 e~l A l [l + 0(e- 2 l A l)j + h e - 2 7l A IA h + 0(e- 2 ^l A l/ 7r ) 

where 7 = 7 (1 — 7/71") -1 and 

1 /vr 2 \ , 7 / 7T 2 
ai = — tan — , b\ = — ~ tan 

7T \27/ 7T \7T — 7 

Thus we formally obtain for 9 — > — 00 

+00 



dt k 
d9 



1-ai ldXe^l^ + dLk 



e e + 0(e 2e ) + 0{e Vf9 l v ) (7.25) 



Whether e or e 7 ' *" is the dominant order depends on who is larger between 
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7 and 7r/3, but it does really affect our argument. In any case self-consistency 
requires that the coefficient of (e ) vanishes, yielding the value of the real part of 
the second integral in eq.(7.23). 

Collecting the results of eqs.(7.24) and (7.25) into eq.(7.22) one obtains in 
explicit form the kink contribution to the energy 

+00 
E P =° ~ - \m 2 (3 cot £ - ^ J d9 cosh 9 [£($) + i(9)] (7.26) 

—00 

The first term is the universal conformal Casimir energy, from which we read the 
correct central charge c = 1 of the mT-sG model. The second term, linear in /3, 
exactly coincides with minus the scaling bulk free energy (see eq.(5.7)). As is well 
known, this is not an accident, since the UV-independent part of the complete 
ground state energy E c + E ~ lim7V->oo En for small r can be calculated through 
Conformal Perturbation Theory (PCT). For the mT model the unperturbed con- 
formal theory is the massless Thirring model, while the perturbation is the mass 
term. This has noninteger dimension for 7 7^ 7r/2, yielding only noninteger powers 
of r in the perturbation series for the free energy, apart from the conformal (3 
term. Hence no term linear in r may appear. This argument has been often used 
to compute the scaling bulk energy of Perturbed Conformal Models from the TBA 
derived from the corresponding scattering theory. In our case we verify the validity 
of the argument since the scaling bulk energy can be directly calculated from the 
microscopic BA solution. 

The last integral in eq.(7.26) represents the resummation of the PCT. The fact 
that it must contain non-integral powers of r, causing non-analyticity at r = 0, 
can be established directly from the original equation (7.18), since this cannot be 
expanded in a Taylor series of r. To establish the absence of any integral power 
of r, and especially of the first power, without referring to PCT, a more elaborate 
treatment of eq.(7.18) is necessary. This shall be the subject of next section. 
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7.4. The NLIE as a Riemann-Hilbert (or Wiener-Hopf) problem. 

In this section we recast the non-linear integral equations (5.9) for the mT-sG 
model as a Riemann-Hilbert (or Wiener-Hopf) problem. The XXZ thermodynam- 
ics can be treated analogously. The techniques involved are rather standard, having 
been used very often to study Bethe Ansatz systems coupled to external fields at 
zero temperature (see for instance [22]). Here we adapt them to the regime T > 
at zero external field. 



Let us define, 



m - hfew*) 



(7.27) 



For large j3, the roots 9^ defined by eq.(5.18) become a continuous distribution 
with a density precisely given by p{9). For arbitrary real values of /3, p{9) is just a 
continuous function that we are interested to find. From eqs.(5.14) and (7.27) we 
find 

dv 



d9 



2np(9) -m/3 cosh9 



(7.28) 



This equation combined with the 9 derivative of eq.(5.21) and with eq.(5.17) yields 
upon partial integration 



p(9) -^cosh9 

2lT 



+oo 



dp Gq(9 — p) 



keTL 



(7.29) 



The asymptotic behaviour given by eq.(5.20) shows that the roots 9k accumulate 

at ±oo. 

The smaller positive root will be called b = 9$. Then, since there are no roots 
in the interval — b < 9 < b, it is convenient to write eq.(7.29) as 

+b 

p(9) - ^ cosh 9- [ dp G (9 - p) p(p) = B(9) (7.30) 

27T J 
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where 

B{6) = f d/i G (9 - ft) 

\u\>b 

For small (3, eq.(5.20) applies and we find 



p(p) -^2$(P- e k 

k£7L 



(7-31) 



/3^0 / 2tt 

= log ( — r 



-oo. 



Hence, for small /3 -B(#) tends to zero and we can consider it a perturbation. 

Eq.(7.30) has the appropriate form to be transformed into a Riemann-Hilbert 

m/3 
~2F 



problem upon Fourier transformation. Notice that the function ?r-cosh# — p{6) 



vanishes for large 9 and can be Fourier transformed as follows: 



+oo 






X+(lu) = e +luh I d6 e~ iue 

+b 


cosh 6 - 

_2vr 


~P{0) 


-b 




X-(u) = e-' iuh J d6e- iue 


cosh 6 - 

_2tt 


-m 



(7.32) 



where the functions X±{u) are analytic for =|=Imc<j > 0. In addition, 



X + (uj) = X_(-w) 



Now, Fourier transforming eq.(7.30) yields 



+b 



mf3 



[1 - G (u)] / e- luM p{e) d6=—- c\i b {uo) + e- tuJb X+(uj) + e +luJb X_(uj) + B(u) 



-b 



(7.33) 



where the function ch^(cj) is defined by 



+1 



ch fe (u;) 



dee~ Ma coshtf 



sinh(l + iu;)6 sinh(l — iu)b 



1 + iu 



+ 



1 — icu 



(7.34) 
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The kernel 1 — Gq(uj) can be factorized as follows 

9 

sinh(^ 



1 - G (u) 



27 I 



27 
1 



sinh(^)+sinh([^-l]7ru; 



K + (lu)K-(lu) 



The function K±{u) is analytic and non-zero for =plm (a;) > 0, K+(u) = K-(—u) 
and explicitly one finds 



K±(u) 



Let us define 



1,1 2 7 



ry\—iu/2 / /y\ 1/2— 4(71-/7— 1)w/2 

tt/ v tt/ r(J-f)ri-(^-i)f 



+6 



/±M 



=fiu)b 



e'~ I d9e- luJV p{6). 



-b 



The function f±(u) is analytic in qTmcJ > 0. Eq.(7.33) yields a pair of equations 
f+ 



K+ 



-K 

-K 



-iub van C h h + X+ e- 2iujb + X. 



■2^cn 6 i-^ + 



e -iub $ 



+ 



,+iujb m/3 
2% 



ch;, + X_ e 



+2icob 



X + + e+^ 6 S 



(7.35) 



Projecting eqs.(7.35) into the half-plane Imu > yields 



K-X- 



K- ( ^ ch & + B) e~ iub 



/v-(^ch 6+J B)e ? 



K_X + e 
X+ X_ e+ 2 ^ 6 



—2iu>b 



0. 



(7.36) 



The projection [. . .]± can be explicitly performed as 

dJ F{uj') 



[Ftp) 



i± 



± 



27ri a; — on 



- , =Flmw > 



Equations analogous to (7.36) follow upon the exchange u <-> —uj, ± <-► =p. 
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More explicitly, the first of eqs.(7.36) reads 



+00 



Kju)Xju) 



AjJ e ~iu'b 

— - -,K(J) 

2iri uj — uj 



^-ch b (u')+B(u') 

In 



+00 



(7.37) 



du' e- 2lJh 



+ / 7T- 7 K ~^') X +^') > Im M > . 

2m uj — uj 



Then, inserting eq.(7.34) in eq.(7.37) and evaluating by residues part of the terms, 
yields 



Kju)Xju) 



im(3 

An 


\jK_{u)-K_{i) 

uj — 1 


e -b K-{u)- 

UJ + % 












+00 

r dw' 

+ J 2m 


-2iu'b 

7 K-(«S) 

uj — uj' 


\x.( J) tmf3 ( 

y ' An \ 


' e b 

-u<- 


- + 
1 


e - b y 




-u' + i) 


—00 










+00 

r dw' 

+ J 2m 


p —iui'b 

7 K-(u')B(u') . 

UJ — UJ 








— 00 



















Setting 



v(u) = Kju) 



Xju) - 



imfi ( e 



An \u — i uj + i 



(7.3* 



(7.39) 



we see that v(u) is analytic in Imw > except for a simple pole at uj — %. We 
then obtain from eqs.(7.38) and (7.39) 



+00 , +00 , 

imf3e b K-(j) f du' e 2iuJ b ,,.,,. f du' e iuj b ,. , - , 
v{u) = — A _ - \, + / — ——, a(u')v(u') + / — : KUu')B(u') 



An u — i J 2ni u + u' 

— 00 



2ni u + u' 



(7.40) 



where 



a[u) = 



K+(u) 



K-{u) 

The ground state energy E can be conveniently expressed in terms of the residue 
of X-{u) at u = —i. In fact, using eqs.(5.8) for rj — ► + , (5.15), (7.28) and (5.17), 
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E may be written as 



E 



+00 
m I dO cosh# 

—00 



P {6) -J2 s ( e ~ °k) 

keTL 



Comparing this with eqs. (7.29) and (7.32), we obtain 



. . u>^— i — %E 
X-(uj) ~ tan 



vr 2 e- b 



mix 27 u + 2 



that is 



£ 



^m/5 h 7r 2 . .. . . 
e cot — hm {uo + i)X-(l)) 

71 27 uj^—i 



We can also express -E in terms of v(—i) 
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E = m cot — 

2 7 



mf3 ine , 



(7-41) 



where 



lim 



tf-M 



7T 



a;^— i Co> 



7^"* r(1 



71 



TV 
2 7 



2 7 



Eqs. (7.40) and (7.41) are exact and provide the starting point to systematically 
compute the high temperature (large b) behaviour. Here we limit ourselves to 
a qualitative analysis: to start we notice that the first term in eq.(7.41) exactly 
provides minus the bulk free energy, which supposedly is the only term regular at 
mf3 = 0. Next we observe that the integration path in the first integral of eq.(7.40) 
can be deformed into the path C+ which runs around the positive imaginary axis. 
Thus we see that this integral gives poles to v(u) in the half-plane Imcu < in 
correspondence with the poles of K + (uj) for Imu > (at u = i, K + (u) has a 
simple zero). These poles are located at Imu = £ n = (27/7?) n, n = 1, 2, . . . , so 
that the integral contributes a series of terms proportional to v(i£ n ) exp(— 2bC, n ) 
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to v(—i) and hence to the energy. Neglecting for the moment the second integral in 
eq.(7.40), we evaluate v(w) iteratively for large b, and conclude that v(—i) has an 
asymptotic expansion in powers of exp(— 26£i) ~ (m/3/27r) 7 ' 7r , including the power 
zero. Inserting this in eq.(7.41) we find that the energy has a term proportional to 
(m/3) -1 plus a series of terms with (m/?) -1 times integer powers of (mfl)^' 71 . The 
effects of the second integral in eq. (7.40) can be taken into account by applying 
the Euler-Maclaurin formula to the definition of B(8), eq.(7.31). To second order 
in the Euler-Maclaurin expansion we then find 



B(u) 



uisinub 

cosuj + ^pW + 



where p(b) can be evaluated to this order by self-consistency This leads to changes 
in the coefficents of the expansion of v(—i) and hence of the energy, but leaves 
unchanged the analytic structure. The same would happen by pushing further 
the Euler-Maclaurin expansion. The structure of this expansion and the powers 
{m(3)^ n '' K agree with Perturbed Conformal Field Theory, since the perturbing 
operator (the mass term in the mT model or the cosine term in the sG model) has 
exactly the scale dimension 2(1 — 7/71). 

8. Comparison against standard TBA in the IR limit m(3 — > 00. 

We can write the non-linear integral equation (5.12) for the massive Thirring 
model as 

Uf = ej — rcosh.9 = —Go* Lf + Gi* Lr (8.1) 

where r = m/3. In Fourier space the two kernels are related by G\(k) = Go(k)e 
so that eq.(8.1) reads there 

i> f = G (-L f + e* k lj) (8.2) 

We have to lowest order 
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• 



L f (9) = log[l + e _e ' W ] ~ e - rcoshe = g{9) 
and thus, to the same first order in the uniformly small g 

O f ~G (-l + e* k )g 

Hence, to second order in g 

L f + L f ~ e e ? - \e~ 2ef + e £fh - \e~ 2 ^ b 

~g(2-u f -u f b)-g 2 (8.3) 

= 2g(l-g)+gW*g 
where the kernel of the convolution operator W has Fourier transform 

2 sinh (£--l)k sinh 2 \h 

W(jfe) = 1 - AG Q (k) sinh 2 ttA;/2 = 1 ±1 ^ 

sinh f |- — g ) fc cosh ^/c 

The traditional TBA is based on the string hypothesis. According to the work 
of Takahashi and Suzuki [2], one needs to write 7/V as continued fraction whose 
structure will fix the type of string configurations allowed as roots of the BAE at 
finite temperature T = 1//3. Standard entropy arguments then lead to a generically 
infinite set of coupled nonlinear integral equations [2]. When 7/7? is a rational 
number, this set of equations can be reduced to a finite system. The simplest 
case is obtained when 7 = n/n, with n an integer larger than 2. One then deals 
with the following system of coupled equations for the real pseudoenergies €j(0), 
j = 1,2, . ..,n 

ei (0) = rcosh9-D*L 2 (9) 

e • = -D* (Lj_i + Lj+i) , j = 2, 3, . . . , n - 3 

(8.4) 
e n _ 2 = £> * (L n _ 3 + 2L„ * 



J n I 



Cra-l = Cn = — -D * Ljj_2 

where Lj = log(l + e~ ej ) and -D(#) = [27rcosh(#)] _1 . The structure of this system 
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is that of the D n Dynkin diagram. In the limit r — ► oo the functions Lj all tend 
to constants, with L\ — > in particular. Since -D(O) = 1/2 we then get the system 
of numerical equations 

e j = -±log(l + e - e - 1 )(l + e-^+ 1 ) , j = 2,3,...,n-3 

e n -2 = ~\ log (1 + e- e - 3 ) (1 + e^-f (8.5) 

e»-i =e n = -ilog(l+e- £ - 2 ) 

with solution 

ej = -log(j 2 -l) , j = l,2,...,n-2 

tn-i = e n = - log(n - 2) 
We now linearize eqs. (8.4) around this r = oo solution, by setting 

Lj = log j 2 + £j , ej = - log(j 2 - 1) - . 2 _ - £j , j = 2,3,...,n-2 ( 8>6 ) 



■2 ■ n - ! 



L n _i = L n = log(n - 1) + 4-1 , e n _i = - log(j - 1) 4-1 

n — 2 



so that the functions 4 fulfill the linear system 



7 2 

' D*(4_i+4+i), j = 2,3,...,n-3 



2 _ 1 ■? _ v J _i ^ ■?" 

4-2 = £>* (4-3 + 2^-1) ( 8 - 7 ) 

n-1 = -D * 4-2 



J 2 -l 

(n-2) 2 
(n- 2) 2 -l 
n — 1 „ 



n-2 



Clearly all 4 are of order g, so that we obtain the order g for the function L\ 
from 

Li ~ e~ €l - \e~ 2ei ~ 2#(1 - g) + gD * £ 2 (8.8) 

since ei ~ r cosh# — log 2 — D * £2- Agreement with the result (8.3) of our TBA, 
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requires 

D*£ 2 = W*g (8.9) 

as soon as 7 = ir/n. This requirement is indeed satisfied: by assuming £2 to be 
given by eq.(8.9) we can easily calculate in Fourier space all other £j's from the 
first n — 3 equations of the system (8.7). The last equation of that system must 
therefore reduce to an identity, and it does. 
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